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Abstract 

This article investigates correlational properties of two-dimensional 
chaotic maps on the unit circle. We give analytical forms of higher-order 
covariances. We derive the characteristic function of their simultaneous 
and lagged ergodic densities. We found that these characteristic functions 
are described by three types of two-dimensional Bessel functions. Higher- 
order covariances between x and y and those between y and y show non- 
positive values. Asymmetric features between cosine and sine functions 
are elucidated. 

1 Introduction 

Knowledge on solvable chaos is useful for designing random number genera- 
tors [U [2J 121 U and Monte Carlo integration [5] . The idea of applying chaos 
theory to randomness has produced important works recently [51 El |H1 12] ■ Geisel 
and Fairen analyzed statistical properties of Chebyshev maps [TU] . They showed 
the mixing properties and higher order moments with higher-order characteris- 
tic functions. Gonzalez and Pino proposed a pseudo random number generator 
based on logistic maps [IT] . Collins et al. [H] have applied the logit transforma- 
tion to the logistic map variable for producing a sequence with a near Gaussian 
distribution. These solvable chaotic properties enable us to design and employ 
chaos for application purposes. 

First, let us consider maps in the form of Chebyshev polynomials of degree 

k 

x t+ i = T k (t t ), (1) 

which map the interval [— 1, 1] onto the same interval. The first few polynomials 
are explicitly T\(x) = x, T2(x) — 2x 2 — 1, and T^(x) = Ax 3 — 3x. Since, there is 
permutability of the Chebyshev polynomials, Tk(Ti(x)) — Tki(x), Eq. (TJJ can 
be expressed as 

x t = T k t(x ). (2) 
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It was shown by Adler and Rivlin that Chebyshev maps with k > 2 are er- 
godic and strongly mixing. This map dynamics has the invariant measure 
u(dx) = — 7= „ . Geisel and Fairen shows that the characteristic function of 

r V ' TT\!1 — X 2 

the Chebyshev maps can be expressed as Bessel function [TU]. They further 
considered the higher-order characteristic function. Following their strategy, we 
consider the characteristic function of two-dimensional solvable chaotic maps 
on a unit circle. We further calculate the higher-order covariance based on the 
characteristic function. 

This article is organized as follows. In Sec. [5J we introduce two-dimensional 
chaotic maps on a unit circle. In Sec. [3j we show that simultaneous covariance 
among two variables is independent. In Sec. [4j we derive an analytical form 
of higher-order covariance among two variables. In Sec. El we compute higher- 
order covariance among two variables with lags. Sec. |6]is devoted to concluding 
remarks. 



2 Two-dimensional solvable chaos 

In this article, we consider two-dimensional maps on a unit circle. Suppose that 
z t = x t + y—lyt denotes a complex number, where Xt is a real number and 
y t is an imaginary part at step t (t = 0, 1, . . .). Then, we define the complex 
dynamics as 

zt+i = z k t , (3) 
where k is an integer. We can also express Eq. ([3]) as 

f %t+i = P k {x u y t ) , 4 - 
1 Vt+i = Qk(xt,yt) ' 

where Pk(x,y) and Qk{x,y) are defined as 

(x + V^ly) k = P k (x,y)+V=lQ k {x,y), (5) 
x 2 +y 2 = 1. (6) 

The first few polynomials are explicitly given by P\{x, y) = x, Qi(x,y) = y, 
P 2 (x,y) = x 2 -y 2 , Q 2 (x,y) = 2xy, P 3 (x,y) =x 3 -3xy 2 , Q 3 (x,y) = 3x 2 y - y 3 , 
Pi{x, y) — x 4 - 6x 2 y 2 + y 4 , Qa{x, y) = 4x 3 y - Axy 3 , Ps(x, y) = x 5 - I0x 3 y 2 + 
5xy 4 , Q 5 (x,y) — 5x 4 y — I0x 2 y 3 + y 5 , P 6 (x,y) = x e — I5x 4 y 2 + I5x 2 y 4 — y e , 
Q 6 (x, y) = 6x 5 y - 20x 3 y 3 + 6xy 5 , P 7 (x, y) = x 7 - 2lx 5 y 2 + 35x 3 y 4 - 7xy 6 , and 
Qr(x, y) = 7x e y - 35x 4 y 3 + 2lx 2 y 5 - x 7 . 

In general, Pk(x, ±Vl — x 2 ) = Tk(x) is satisfied. Specifically, Qk{x,y) for 
odd ordered k is equivalent to Qfc(±yl — y 2 , y) = — Tfc(y). 

If we set an initial condition zq = xq + >/—Tyo on the unit circle \zq\ — 1 , Zt 
is also mapped on the unit circle. In this case, Eq. ([5]) can be rewritten as 

exp(V^ie) k = exp(fc(9V^T), (7) 
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where 9 denotes the argument of (x, y) on the two-dimensional plane. It is 
convenient to represent the polynomial Pk(x,y) and Qk{x,y) in the form 




Figure 1: 800 steps of a trajectory of the two-dimensional chaotic map for k = 2. 
The initial value is given by (x , y ) = (-0.820000, 0.572364). 

By introducing 9 t as the argument of Zt, we have 

9 t+1 = k9 t . (9) 

The solution of Eq. ([9]) can be written as 

9 t = k%, (10) 

by using 6q, denoted as the argument of zo- Therefore, Zt = Xt + \[— Tyt is 
rewritten as 

z t = cos(/c 4 6l ) + ^/^lsinik^o) = exp(jfe*^ V = T). (H) 

Eq. (|T0)) is ergodic and has the constant invariant density pe (<?) — (0 < 
9 < 2ir) since Eq. ((9]) is a Bernoulli map on mod 2ir. 

Transforming the orthogonal coordinate (x, y) into the polar coordinate (r, 9) 
by x = rcos# and y = rsinfl, we have pn(r) = 8(r — 1). Therefore, the joint 
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invariant density of x and y can be described as 

/ \ lns , v 0(9, r) S(Jx 2 +y 2 -l) 

where S(-) represents Dirac's S- function. The marginal density in terms of x is 
given by 

Px{x) = \ p X Y(x,y)dy 



dy 



-l 

i r 1 



27T J_i yjx 2 + y 2 
1 z-^+i- 1 ^(f) 



I / 

7f i|x|-l 



|x|-l v/(t+l) 2 -x2 
1 



di 



71- yr 



In the same way we obtain 



1 1 



Pv(y)= pxv{x 1 y)dx= — (13) 
J-i 7r^l-y 2 

Note that px{x) and py{y) are the same as the ergodic density of the Chebyshev 
maps. 



3 Simultaneous covariance 

Next, let us consider auto-correlations of x and y and cross-correlation between 
x and y. Obviously, mean values of x and y are given as zero. 



1 T_1 Z" 1 Z" 1 x 
l im = / zpx(a;)da;= / da; = 0, (14) 



t=o 

T-1 r l -l 



lim = / / — 7 ===dy = 0. (15) 

We shall introduce four types of correlations: 

cxx(t) = lim — V* x t x t+T = / da;/ dya; P fe o • • • o P k (x, y)p X Y (x, y) 



T-1 



(16) 



cyy(t) = lim ^ V y t y t+r = / da; / dyyQ k o---oQ k (x,y)p X Y(x,y) 
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(17) 

cxy(t) = lim — \^ x t yt+T= I dx dyxQ k o • • • o Q k {x,y)p XY {x,y) 
T ^ T tZ J-i J-i ' v ' 

(18) 

cyx(t) = lim^ — 22y t x t+T = dx J dyy P k o ■ ■ ■ o P k {x, y)pxy(x, y) 

(19) 

Transforming the orthogonal coordinate (x, y) into the polar coordinate 
(r, 9), we can calculate Eqs. (fl"6|) to (Tl9t as 

cxx(r) = / cos0cosfc T 0d0 = i* ltfcT (20) 

2tt Jo 2 

1 /" 27r 1 
cyy(t) = — I sin9shik T 9d9 = -S lk , (21) 
2?r Jo 2 



c XY (r) = — / cos sin fc T 0d0 = (22) 

27T 



i r n 

cyx{t) = — \ sin0cosfc T 0d0 = (23) 
27r Jo 

These are extensions of Chebyshev maps derived by Geisel and Fairen to the 
two-dimensional map |10j . Therefore, the auto-correlations of x and y decay 
for t > 1, and the cross-correlations between x and y are zero. Furthermore, the 
correlation between z t and z t + Tl where 7 is denoted as the complex conjugate of 
•, is also zero, 



lim h z t z t+r = c X x(t) - cyy(r) + V^1(cxy{t) + c Y x(t)) = 0. (24) 

t=0 

Note that Eqs. (|2"0")) to (|2"3"|) are derived by means of the permutability of z fe 
and the orthogonality between P k (x,y) and Q k (x,y). Clearly, from Eq. ([3]) we 
can prove the permutability of z k such as (z k ) 1 = z kl . For k > 1 and Z > 1, we 
also have the orthogonal relations among P k (x,y) and Q k (x,y) 

J dx J dyP k (x,y)Pi(x,y)pxY(x,y) = J cos(Zc0) coa(W)d9 = ^S k 



1 f 27r 1 



(25) 

dx I dyQ k (x,y)Qi(x,y)pxY(x,y) = — I sm(k9)sm(W)d8 = -Sk,i, 

-1 J-l 27T J 2 

(26) 

J dx J dyQ k {x,y)Pi (x,y)p X Y(x,y) = / sin(fc0) cos(Z0)d0 = 

(27) 
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4 Simultaneous higher order covariance 

Let us consider the characteristic function of the simultaneous joint density 
Pxy(x,u), defined as 

1 T_1 

t=0 

/oo poo 
/ e^ ux+v ^p XY (x,y)dxdy. (28) 
-oo J — OO 

Inserting Eq. flT2J| into Eq. (|2"g]l. we have 



OO 

OO 



dO I rdre^ 

2nr 



l(ti cos sin 8) jj 

27rr 

_/ e v^T(«cose +t -sine) d6) = (29) 
27r Jo 

where JP' q (u,v) is defined as 

1 raw 

JP'"(U,V) = — e V^(uco S (p0)+vsin( q e)-n0) dg _ ^ 

27r Jo 

This is similar to the two-dimensional Bessel function which was studied by 
Korsch et al. [T3], however, it is a bit different from it. They define the two- 
dimensional Bessel functions with three integer indices n, p, and q as 

>'«(«,«) = —/ e V=T(usin(pfl)+usin(9e)-nfl) d g ( 31 ) 

27r J-7T 

In his definition, the two-dimensional Bessel function consists of two sine func- 
tions. However, in our definition this consists of cosine and sine functions. 
Clearly, both the two-dimensional Bessel functions satisfy 

J o M (u,0) = Jo(u), J*>\0,v) = J (v), (32) 

%'\u,0) = J (u), J M (0,u) = J (u), (33) 

where J n {u) is the Bessel function defined as 

J n («) = — / e V^TM-«sin0) d6) _ (34) 

In the one-dimensional case, Eq. (I2U1) is equivalent to the characteristic 
function of Chebyshev polynomials, which is derived by Geisel and Fairen [TU] . 



G 



We can further expand <fr(u, v) in terms of u and v, 
$(u,u) = — V 1 , / (wcos6» + wsin6»)"d6i 



n=0 



1 " (^ZT)n " 

2tt ^ n! ^ 

n=U m— 



/ n 








V m , 





2tt 

cos m 0sin n - m 0d0. 



Therefore, we have 



T-l 



do: / &yx m y n ~ m pxY(x,y) 



1 

= — cos m <9sin"- m (9d6>. (0 < m < n). (35) 
2tt Jo 

We also have the equality 
where B(a, b) denotes the beta function, defined as 



(36) 



1 

B(a,b) = I T^il-rf^dT, (37) 
Jo 

and r(a) represents the gamma function, defined as 

POO 

T(a) = / e~ T r Q ~ 1 dr. (38) 



Inserting Eq. (|36|) into p = m/2 + 1/2 and q = (n — m)/2 + 1/2 and using 
symmetry of cosine and sine functions and T(n + 1) = n\, we obtain 

(- 2r(^±i)r(^) _ (m-l)!!(n-m-l)!! / . x 

|jrayii-raj _ J 27rr(f +1) ~~ n<! V 1 ' " l ' eveIi J _ (39) 

[ (otherwise) 
Hence, the characteristic function of pxy{x,y) is described as 

This is a natural extension of the Bessel function of degree to the two- 
dimensional case, 

*<•>-£ (5^- < 41 > 
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Since we can further calculate the m-th order moment of Xt and the n — m-th 
order moment of yt as 

(X m ) = J Ax J dyx m p XY {x,y) 

1 /-27T ( (m-1)!! / \ 

1 / main — rr - (m : even) 

= — I cos'" 6AQ — { n mU ) J (42) 

2vr J \ [m : odd) v ; 



and 



, V'-" , ' I I dx [ dyy n - m p XY (x,y) 



-1 



1 p2ir \ (n— m— 1)!! / \ 

= J_ / sin»-™M0= («-^ :even ) , (43 ) 

27r J I {n — m: odd) 

where m!! = 2 • 4 • 6 • • • m for even m and m\\ = 1 • 3 • 5 • • • m for odd m, we get 
Cov[J^ m Y n ~ m ] (x m y r?l_m ] (J^ m ) (l^ n ~ m ) 



(m-l)!!(n-ro-l)!! 



1 



(),!! 



m!!(n— m")!! 



(m,n : even),^. 
I (otherwise) 

Here, we consider the negativity of even ordered moments. Hammersley sug- 
gested that antithetic variables are effective for variance reduction in Monte 
Carlo integrations [14] . The antithetic- variates method permits estimates through 
the use of negative correlated random variables faster than independent random 
variables. Let us confirm the sign of Eq. ([44|). We get 

1- .., ... = 1 [2> ' 



m\\{n-m)\\ (f)!(^)! 



= 1 - ( | ) < 0, (45) 



since from the definition of combination, we have 

4 ) = (m)I(2=ni)I ^ L ^ 

The equality is satisfied if and only if m = or m = n. Note that Eq. (|44| is 
independent of a value of k. 

Therefore, Eq. (|44p implies that xt and yt do not have any correlations for 
the odd-ordered moments, however, do have a negative covariance for the even- 
ordered moments. Fig. [2]shows the relationship between n and Cov[X m , Y n ~ m ] . 
It is confirmed that the covariance monotonically increases and approaches to 
zero as n increasing. 
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Furthermore, we calculate covariance between x™ and x" 171 , and between 
y t m and y t "~ m - From Eqs. g2|) and (g3]), we have 



Cov[x m ,x n - m ] = Cov[r m ,y 



m \ / n — m 

n—m 



> (n, m : even) 
(otherwise) 

(47) 



The non-negativity of Eq. (|47|) is proven as follows. Let us consider the case 
that n is even. From 

(48) 
(49) 



[l + x) n = {(1+x)*}' 



one has 



E 



n 
m 



x 



2 



2 



m— m— 

Comparing a; m 's coefficient, we get the following inequality 

2 



n 
m 



> 2 

2 



Therefore, we obtain 



f 

2™ 



m \ / n—m 

n—m 

1 



2" 



??? 



(50) 



-( m / 




/ 11 




)-( 


m 













> 0. 



(51) 



5 Higher order covariance with lags 

More generally, we can introduce a characteristic function of the joint density 

— n 
1-8 

T-l 



between and Vt+q 1 



V XY (u,v) = lim - V e V=i(ux t+p +v Vt+ ,) 

t=0 

(exp(y^T(u Pfc o ■ ■ ■ o Pfcj z, y) + f Qfc o ■ ■ ■ o Qk jx, y))) ^ 



i ,i 

dx / dyexp(v / ^T(wft ° ■ ■ ■ o P k (x,y) + v Q k o • • • o Q fe (a;, y)))pxr( 



1 

2^ 



l ./-l 

2-n 



P q 
k p ,k" , 
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Figure 2: The relationship between n and Cov[X m , Y n m ] for k = 2. 



Similarly to <&(u,v), from the expansion in terms of u and v, we obtain 



1 T_1 r 1 f 1 



k o • • • o Q k (x,y)p X Y(x,y) 



1 

2Wo 



cos m (fc p 6>) sin n - m (fc«0)d0. 



(53) 



By using 



cos m (F6>) sin n - m (Jfe 9 0) 

2 m (2 N / = T)«- 



Ont./ i~\ri— m ^— < ^— < ^ ^ 



n — m— s 



( n m )' ^l\(2r-m)k p + (2s-n+m)k q 



r!(m — r)! s!(n — m — s)! 



(54) 



and 



we obtain 



1 



(55) 



T-l 
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c-ir 



E™ y^n-m m ; (n—m)\_ /_-. \_ s r 
r=0 Z^s=0 H^m^Ty! s!(n-m-s)! V ^ o 0,(2r-m)fcP + (2s-n+m)fc9 



2" 







(m, n : even) ■ 
(otherwise) 

(56) 



Since we further have 



dx 



1 

2Wo 



Pfe o ••• o P k {x,y) 



Pxy(x,v) 



cos m {k p 9)d6 



and 



(m-1)!! / \ 

m ii («* : even) 
(to : odd) ' 



(Y t n + q m ) = I dx I d y Qk ° • • • ° Qfc(*. y) w) 



i /- 2,r 

— / cos"- m (fc«6i)d6i 
2tt 7o 



(57) 



'0 

(n— m— 1)!! / \ 

—, TTf— (n — m : even) 

(n— m)!! V / 

(n — m : odd) 



(58) 



we get 



Cov[x™ p ,y t " + - m ] = <* t + P ^ t +~ m > - (^ t + P )(y t +7 n ) 

(— 1) 2 y^^i Y^n— r, 
2™ 2^r=0 2^is=0 



(m-l)!!(n-m-l)!! 
m!!(n — m)!! 



m 
r 



n — to 
s 



( — 1) S ^0,(2r-m)fcP + (2s-n+m)fc9 







Kohda et al. showed that the higher-order covariance of Chebyshev maps have 
no correlation [TS]. We use their derivation in our case. According to Kac's 
statistical independence |16j when in Eq. ([5£ 



(2r-m)k p + (2s-n + m)k g = 0, (0 < r < to;0 < s < n - to) (60) 

holds for any k p and k q if and only if r = to/2 and s — (n — to)/2, k p and /c 9 
are called linearly independent. Then x™ +v and yt+™ are statistically indepen- 
dent [T5] , 

Let consider the case that to and n are even. From elementary facts about 
the theory of numbers, we know that 



(to, n : even) 
(otherwise) 

(59) 



N = k e + r (0<r<k), 



(61) 



11 



where TV is a natural number, and k, e and r are non-negative integers. In the 
case that 2r — m > 0, 2s — n + to < 0, and p < q we have 



(2r-m)k p + (2s-n + m)k q = {(2r - to) + (2s - n + m)k q - p }k p 

= {(k ei +r')-(k e2 +s')k q - p }k p 
= (k ei +r'- k E2+q - p - s'k q ~ p )k p . (62) 

Therefore, if [(2r — m)/fc] =0, [(2s — n + m)/fc] = 0, and e\ = e2 + q — p hold 
then (2r — m)fc p + (2s — n + m)k q = is satisfied for integers r and s other 
than r = m/2 and s = (n — m)/2. When to > k, and n — to > fc, we have 
\{2r — to) /fc] = and [(2s — n + to) /k] = 0. Therefore, m > k ei and n — m> k e2 
would be satisfied. Namely, when n < k ei + k e2 = k e2 (k q ~ p + 1), x"^_ p and 
yf+™ are statistically independent. This implies that q — p goes infinity, x™ +p 
and y"_)T g m become statistically independent in an exponential manner. 

Fig. El shows Cov[X% p) Y t n + - q m ] for (p,q) = (0,1), (0,2), (0,3), (0,4), (0,5), 
and (0, 6). As shown in figures, we found that the covariances decrease as \p — q\ 
increasing. The range of the covariances approach to zero as q increasing. 

Obviously, Eq. (|60|) has solutions r = m/2 and s = (n — m)/2. A sum of 
the contributions for r = m/2 and s = (n — m)/2 in Eq. (fTTj) is equivalent to 
^ m ~m\\(n-m)\^ A 1 Since Cov[A t ™ p , Yj""" 1 ] is less than zero from the numerical 
simulation, for solutions other than r = m/2 and s = (n — m)/2 of Eq. ([6H)1 . 
it should satisfy that a sum of negative contributions is greater than a sum of 
positive contributions. 

We may consider two types of second-order characteristic functions with lags. 
Note that Geisel and Fairen [10] considered a similar second-order characteristic 
function for the Chebyshev maps. Their characteristic function corresponds to 
fxx(i,«) in our definition. 

T-l 

l(uxt+p+vx t + q ) 



^xx(u,v) = lim ^ V , 

T— ¥00 1 £ — ' 

t=0 

= _L I" 2 * e V=T(ucos(k*>8)+vco S (ki8)) d Q (gg) 

2tt Jo 

T-l 

^Yy(u,v) = lim - V e V=l(uy t+p +vy t+q ) 
t=0 

2?r Jo 

Similarly to \P xy(u, v), from the expansion in terms of u and v, we obtain 

*xx(u,v) = jrt±)lf2 ( " ) (XP +p X^ q m )u m v n ~ m , (65) 



n— m— 



* YY (u,v) = ]T L^L J2 ( n m ) {Y™ p Y^™)u m v n - m , (66) 



n=0 m=0 
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Figure 3: Scatter plots of Cov[X t ™ p , F t ™~ m ] in terms of n (0 < m < n) at 
k = 2 and p = 0, (a) q = 1, (b) q = 2, (c) q = 3, (d) q = 4, (e) q = 5, and 
(f) q = 6. Filled squares represent theoretical values, and filled circles values 
obtained from numerical integration. 



T-l 



where 

\ t+p t+q I T ^ oo T £^ 
t=Q 
T-l 

\ I t+p 1 y+q I - rp Vt+pVt+q 



m n-m 
^t+p^t+q 



1 

2^ 



cos m (k p 6) cos n - m (k q 6)d8, 



1 f 2 * 

— / sm m (k p 9)sm n - m (k c '9)de. 
2tt J 
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By using 

cos m (fc p 6>) cos n - rn (k q 6) 

1 / V^Tfc p 8 , c -v^Tfc p eyn 1 / s/=lk q 6 , -V^Ifc^n-m 

2m / 2 n— m 

m n—m . / 

J_ w! l n ~ m J ! V=T[(2r-m)fc p + (2s-n+m)fc 1 '10 

2« r i( TO _ r )i s u n _ TO _ s \\ 

r=0 s=0 V ; V ' 

sin m (P(9) sm n - m (k q 9) 

1 / sf^\k v 6 _ e -sf^lk v d\m \ _ e -V=Ifc 9 0\n-m 



2v^T) m (2v^T)™- 

I ,n m n—m . , •.. 

(-L) 2 x ^ m< (71-771)1 ^-r-s c ^=l\(2r-rn)k' , + (2s-n+m)k q 

2" 2-^ t\{jtl — r)! s\(ti — m — 5)! ' 

r=0 s=0 \ I \ I 



J_ m n-m 



therefore, we have 

{1 v-^ m sr^n—m m' (n—m)! r / \ 

F Lr=0 L s =0 r!(m-r)! s\(n-m-s)\ 0,(2r-m)kP + (2s-n+m)ki (TO, n : even) ^ 
(otherwise) 

{(— 1) "5" y-^ m v-^n— m m' (n-m)! / -,\_ r _ S r / \ 

2" 2^r=oZ^=0 r!(m-r)! s!(n-m-s) ! ^ °0,(2r-m)kP +{2s-n+m)k« K™>, " : even J (68) 

(otherwise) 

We further have 

U0V[A t+p , A t+? J - {A t+p Ji t+q ) - {JL t+p ){Ji t+g ) 



1 ^ m \^n— m I 771 \ I 71 ffl i ^ 
W l^r=0 l^s=Q \ r j \ S ' °0,(2r-m)kP + {2s-n+m)ki 



(69) 



/(m-l)!!\ / If 

( - m!! j (to, n: even) 

(otherwise) 

(70) 

A sum of contributions for r = to/2 and s — (n — m)/2 in Eq. ([TO")) is 
equivalent to ( ^ m m! |- >!! ) 2 . If Eq. (1TO1) has other solutions than r = m/2 and 
s = (ra — to)/2, then the covariance positively increases. Therefore, we could 
prove Cov[A-™ p ,Xr + - m ] >0. 
We also have 

Cov[r™ p ,r/;-™] = (r™ p r t r 9 m > - TOWD 

-l)f y^m ^n-m / 771 \ ( 71 — 771 \ , x_ r _ s c 

~ Lr=0 Ls=0 V T J I S )\~ L ) <J0,(2r-m)kP + (2s-n+m)ki 

2 



(n-m-l)!!\ / W 

(n-m)!i j (to, n: even) 

(otherwise) 



(71) 
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Fig. [5] shows covariance between X^_ p and X™ + ™, and between Y™ p and 
Y^T q m . It is found that Cov[X£+ , X™+™] shows non-negative values, and that 
Cov[Y™ p ,Y t \- q m ] shows non-positive values. We found that Cov[X™ p , F 4 "~ m ] 
takes the same non-positive value as Cov[Yj" l p , Y^"" 1 ] for p ^ q from Figs. [3] 
and[5] The reason is because cos m (fcP(9) sin n - m (> 9 0) and sin m (fcP(9) sm n ~ m (kW) 
have the same area to the x-axis, but cos m (k p 8) cos n ~ m (k q 9) is different from 
them as shown in Fig. |4j 

A sum of the contributions for r — m/2 and s = (n — m)/2 in Eq. (|71[) is 
equivalent to ( ) 2 . Since Cov[F t " l p . Y^ q m \ is less than zero from the 
numerical simulation, for solutions other than r = m/2 and s = (n — m)/2 of 
Eq. (|60p . it should satisfy that a sum of negative contributions is greater than 
a sum of positive contributions. 




Figure 4: The wave forms of cos m {kP8)sin n ~ m (kW), sin m (fcP<9) sm n - m (kW), 
and cos m (fc p 6>) cos"- m (fc«6>) for p = 0, q = 1, n = 10, and m = 4. 

Therefore, it is suggested that ^xx(u,v) ^ fyy(u,i)) ^ %y(H,«) for 
q =/= p from numerical simulation. This also implies that three types of two- 
dimensional Bessel functions are not equivalent; 

i /■2tt 

JP> q (u,v) = — / e V^T(«cos(p9) +t ,cos( 9 e)) d6/; (72) 
27r Jo 

1 /-2w 

J™(u,v) = — / e V=T(«sm(>0)+i;co S ( g e)) d6 , ) (73) 

27r Jo 

1 /-27T 

JS«(U,«) = — / e V=l(u S in(p9)+vsin( q 9)) dd _ (j^ 

2tt Jo 
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Figure 5: Scatter plots of Cov[X™ p , X?+™] and Cov[F t ^ p , Y£ ~ m ] in terms of n 
at k = 2 and p = 0, (a) q = 1, (b) q = 2, (c) q = 3,(d)q = 4, (e) q = 5, and (f) 
q = 6. Unfilled squares represent theoretical values of Cov[X]^ p , X£^ m ], filled 
squares numerical values of Cov[X™ p , X™+™], unfilled circles theoretical values 
of Covfr™ p , Y^™}, and filled circles numerical values of Cov[Y™ p , Y^T p m ]. 



(f) 



6 Conclusion 

We studied two-dimensional chaotic maps on the unit circle, which is an ex- 
tension of the Chebyshev maps to two-dimensional map on the unit circle. We 
examined correlational properties of this two-dimensional chaotic map. We gave 
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analytical forms of higher-order moments. Furthermore, we derived the charac- 
teristic function of both simultaneous and lagged ergodic densities. We found 
that these characteristic functions are given by three types of two-dimensional 
Bessel functions. We proved four theorems and proposed two conjectures as 
follows: 

Theorems: 

1. The higher-order covariances between xt and yt shows non-positive values 
for integers n and m (0 < m < n): 

Cov[X m ,Y n - rn ] < 0. (75) 

2. The higher-order covariance between Xt and Xt shows non-negative values 
for integer n and m (0 < m < n): 

Cov[X m ,X n - rn } > 0. (76) 

3. The higher-order covariance between y t and yt shows non-negative values 
for n and m (0 < m < n): 

Cov[Y m ,Y n - rn ] > 0. (77) 

4. The higher-order covariance between Xt+ P and Xt+ q (p ^ q) shows non- 
negative values for integer n and m (0 < m < n): 

Cov[X^X^ g m }>0. (78) 

Conjectures: 

1. The higher-order covariances between Xt+ P and yt+q (p ^ q) shows non- 
positive values for integers n and m (0 < m < n): 

Cov[xr +p ,Y t n - q m }<0. (79) 

2. The higher-order covariance between yt+ P and yt+q (p ^ q) shows non- 
positive values for n and m (0 < m < n): 

Gov[Y t - ,Y t l~ m }<0- (80) 

Therefore, we can generate antithetic sequences as xq, yo, x%, y±, . . . , xt, yt, ■ ■ ■ 
or yo,yii 2/2) • • • j Dti ■ ■ ■ obtained from Eq. Asymmetric features between 

cosine and sine functions were elucidated. Using the proposed two-dimensional 
chaotic map, we can generate antithetic pseudo random sequences for Monte 
Carlo integration. 
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